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ABSTRACT 


A preliminary investigation was conducted to evaluate the lateral slosh 
characteristics of a 23-cm- (9-in. -) diam mercury propellant tank with a 
positive-expulsion diaphragm and 17. 5% ullage. Data are presented for tank 
sinusoidal acceleration levels between 0. 05 and 0. 5 g at frequencies ranging 
from 5 to 50 Hz. Results indicate that the slosh characteristics are highly 
nonlinear at acceleration levels approaching those in flight, and depend 
heavily on the shape and stiffness of the diaphragm. Nyquist plots of driving 
force over acceleration are shown to be a useful tool for determining natural 
frequency, damping, and modal mass characteristics. A computerized non- 
linear least- squares method for extracting the modal parameters from the 
Nyquist plots is described and results of applying the method are presented. 
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I. INTRODUCTION 


The Solar Electric Propulsion (SEP) Thrust Subsystem Engineering 
Model Technology Program is aimed at the design, development, and func- 
tional test of an engineering model of a solar electric propulsion subsystem 
for a SEP spacecraft. A portion of the Program is centered on the tank that 
contains the mercury propellant for the ion engines. Current emphasis is on 
a spherical tank equipped with an elastomeric diaphragm for the positive 
expulsion of the mercury. 

An important part of the development of the tank design is the evaluation 
of the slosh characteristics of a partially off-loaded tank. This report 
describes the results of a study of the lateral slosh characteristics of a 
23 -cm- (9 -in. -) diam prototype mercury propellant tank developed during the 
Solar Electric Propulsion System Technology (SEPST) Program (Ref. 1). 

The general configuration of the tank is indicated in Fig. 1. 

II. TEST OBJECTIVE AND SCOPE 

The primary objective of the test program was to gain a preliminary 
understanding of the dynamic slosh characteristics of the prototype tank. 

Of particular interest was the dynamic behavior of the fluid at vibration 
levels approaching those anticipated during the launch of an SEP spacecraft 
aboard a Titan-IIID or shuttle launch vehicle. A secondary objective, estab- 
lished during the course of the study, was the development of an improved 
method for reducing the experimental slosh dynamics data and constructing 
equivalent mechanical models of the sloshing fluid. 

Early in the study a number of items were identified as important 
slosh parameters. These are listed below: 
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(1) Gravity field (1 _g vs quasi-static launch). 

(2) Tank diameter. 

(3) Ullage volume. 

(4) Ullage fluid (gas vs liquid). 

(5) Ullage pressure level. 

(6) Diaphragm stiffness (thickness and material). 

(7) Dynamic excitation level. 

(8) Dynamic excitation direction (vertical vs lateral). 

Because of programmatic constraints it was necessary to limit the 
study to a single tank configuration and ullage volume. Emphasis was placed 
on establishing the general slosh behavior as a function of excitation level, 
and, in particular, evaluating the natural frequencies, damping, and modal 
masses of the lateral slosh modes at sinusoidal excitation levels between 
0. 05 and 0. 5 j* peak. It was assumed that follow-on work would be performed 
to evaluate the slosh problem further. Determining the scope of future test 
requirements and the generation of general computerized data reduction tech- 
niques compatible with future testing requirements was considered an impor- 
tant part of the current study. 

The test configuration is summarized below: 

(1) Tank: 23-cm-diam sphere. 

(2) Diaphragm: 1 . 6-mm-thick neoprene. 

(3) Fluid: mercury. 

(4) Ullage volume: 17. 5%. 

(5) Ullage fluid: air. 

5 2 

(6) Ullage pressure: 10 N/m (1 atm). 

III. APPROACH 

The problem of determining the dynamic properties of a sloshing fluid 
is basically one of system identification from measured response data. The 
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usual approaches are based on the measured amplitude response of the fluid 
to steady- state sinusoidal base excitation. 

Using a technique described in Ref. 2, the sloshing fluid was repre- 
sented as a rigid mass plus a number of uncoupled single-degr ee-of-fr eedom 
oscillators, each representing a lateral slosh mode. Each oscillator is 
described in terms of its natural frequency, damping ratio, and mass. The 
problem of determining the dynamic properties of the sloshing fluid was thus 
reduced to determining the natural frequency, damping, and mass associated 
with each dominant slosh mode. The rigid "residual" mass of the fluid 
equals the total fluid mass minus the sum of the slosh masses. 

Though methods for determining the dynamic properties (natural fre- 
quency, damping, and mass) of multiple -degr ee-of-fr eedom systems are 
generally based on system amplitude response to steady-state sinusoidal 
excitation, it is desirable also to make use of the phase relationship between 
the driving force and the system response. When both amplitude and phase 
data are available, it is possible to generate Co-Quad or Nyquist plots that 
greatly aid the determination of the dynamic properties. 

In the current program the propellant tank was excited with a constant 
amplitude sinusoidal base acceleration at frequencies ranging from 2 to 50 Hz. 
The driving force at the base, the base acceleration, and the phase between 
them were recorded. These data were plotted in the form of Nyquist phase- 
plane plots. The dynamic properties were then extracted from the Nyquist 
plots using computerized curve fitting techniques. The data reduction tech- 
nique is described in greater detail in Section VI. 

IV. APPARATUS AND INSTRUMENTATION 
A. Propellant Tank 

The tank used in the slosh test was a modified version of the 23-cm- 
diam SEPST prototype tank described in Ref. 1. To allow the slosh motion 
to be observed, the tank was fitted with a clear Lucite upper hemisphere and 
a grid pattern was painted on the neoprene diaphragm. 
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B. Sloshing Apparatus and Instrumentation 

To provide for near-frictionless lateral motion, the tank was supported 
on a carriage free to roll on linear ball bushings. The tank/carriage assem- 
bly was driven by an Unholtz-Dickie Model 4 electrodynamic shaker with a 
10-cm (4-in. ) stroke and 670-N (150-lb) force capability. An accelerometer 
was positioned on the carriage for acceleration measurement, and a precision 
load cell was positioned between the shaker and the carriage for force mea- 
surement. The initial configuration with two shakers is shown in Fig. 2. 

Only one shaker was used during most of the testing. 

Signals from the accelerometer and load cell were fed into the shaker 
control and signal processing console shown in Fig. 3. Both signals were 
conditioned through a tracking filter to strip off response harmonics not at 
the driving frequency. The amplitude and phase data were then printed on 
the teletype along with the driving frequency. 

All moving components of the shaking apparatus were carefully weighed. 
Calibration on the entire system was then verified by initially shaking the 
tank in a completely filled condition, and noting that the force over accelera- 
tion agreed with the measured system weight, and that the phase angle was 
zero. 


V. TEST RESULTS AND DISCUSSION 

A. Preliminary Test With Water 

Preceding the system verification test with the tank completely filled 
with mercury, a test was conducted using water in place of mercury. This 
test was primarily aimed at debugging the test procedures and equipment 
before introducing the certain degree of hazard associated with the mercury. 
However, the test using water also provided some significant observations. 

An 82. 5% filled (17. 5% ullage) condition was used in the water slosh 
test. This is the same level used later in the mercury tests. As shown in 
Fig. 4, the diaphragm was generally spherical in shape and did not have the 
concave dimple exhibited when the water was replaced with the same volume 
of mercury. The shape of the diaphragm when 82. 5% filled with mercury is 
essentially that shown in Fig. 5 for the 20% ullage condition. 
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In both cases, water and mercury, it is clear that the diaphragm sig- 
nificantly alters the free surface of the fluid by supporting the fluid away 
from the tank wall. As a result the classical first free surface slosh mode, 
which is typified by a general fluid pendulum motion, does not occur. At low 
vibration levels (less than 0. 1 gj the system stiffness appears to be almost 
entirely controlled by the stiffness and shape of the diaphragm. The first 
mode shapes for the water and mercury were entirely different. The first 
slosh mode for the water was a first cantilever mode motion of the entire 
free surface hemisphere at approximately 27 Hz, over ten times the theo- 
retical resonant frequency for the same fluid volume with a free surface. 

The first slosh mode with the mercury was a wave action mode asso- 
ciated with a rolling of the diaphragm dimple from side to side at about 
2. 5 Hz. However, the dominant mode was the second mode, which corres- 
ponded to flexing of the side walls of the diaphragm at about 9 Hz. The 
mercury slosh results are described in greater detail later. 

From a comparison of the slosh behavior of the two fluids with the 
theoretical free surface behavior, several conclusions can be drawn. 

(1) Theoretical models available for free surface sloshing do not 
appear to provide a reasonable means of predicting slosh modes 
and frequencies with the tested diaphragm. 

(2) Because the slosh modes and frequencies are apparently highly 
dependent on the diaphragm stiffness and deformed shape, the 
slosh frequencies should not depend on the vertical gravity field 
in the classical Jg manner. However, the difference between 
the water and mercury mode shapes and frequencies suggests 
that the slosh behavior could be a strong function of the acceler- 
ation field strength. 

(3) Because of the inappropriateness of the free surface slosh theory, 
direct means for extrapolating to the behavior of different ullages, 
larger tanks, different diaphragms, and different acceleration 
fields are not currently available. As a result emphasis will 
have to be placed on acquiring empirical slosh data. 
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B. 


Mercury Slosh. Test (17. 5% Ullage) 


Following the system verification tests, the tank was filled with 67 kg 
(147 lb) of mercury to give an 82. 5% filled condition. Tests were run using 
constant acceleration sinusoidal base excitations of 0. 05, 0.1, 0. 3 and 0. 5 g 
peak, at frequencies ranging from 2 to 50 Hz. Though the slosh was only 
moderately sensitive to acceleration level at both acceleration extremes, 
there was an abrupt change in the system behavior between 0. 1 and 0. 3 g_. 

This change is indicated in the quadrature response plots presented in Fig. 6. 

The quadrature response plots are plots of the normalized component 
of the driving force, which is 90 deg out of phase with the driving acceleration. 
Since the force associated with driving the rigid tank and nonre sonant fluid 
mass is zero or 180 deg out of phase with the driving acceleration, it does 
not influence the quadrature response. For this reason the peaks in the 
quadrature response more accurately represent the natural slosh frequencies, 
than do peaks in the total response. In Fig. 6 the quadrature force is nor- 
malized by dividing by the total fluid mass times the acceleration level. The 
normalized quadrature force can thus be considered as a force amplification 
factor relative to the reaction force required if the fluid were a rigid mass. 

The mode shapes observed to correspond to the resonant peaks in 
Fig. 6 have been sketched in Fig. 7. Because of the inherent uncertainties 
in visually observing mode shapes, the sketches are intended to provide only 
a qualitative insight into the slosh phenomena. As indicated, the first slosh 
mode at the lower vibration levels differed significantly from the first mode 
at the higher levels. At the 0. 05- and 0. 1 -jg levels the first mode consisted 
of a surface wave associated with a rolling of the diaphragm dimple from 
side to side. This mode was highly unstable and quickly coupled with liquid 
swirl. At acceleration levels above 0. 2 g the diaphragm was unable to 
maintain its dimpled shape as the frequency was raised above 3 Hz. At this 
frequency the tank entered a completely new first mode resonance associated 
with a deformed diaphragm shape. In this mode a large fraction of the con- 
tained mercury '’flopped'’ from side to side, stretching the diaphragm and 
impacting the wall of the tank. The mode did not exhibit liquid swirl, but 
became unstable above 5 Hz and disappeared with a jump back to the original 
dimpled diaphragm equilibrium shape. 
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Above 5 Hz the slosh response was similar at all excitation levels with 
a general trend toward higher damping and lower natural frequencies with 
increasing excitation level. As indicated in Fig. 7 the higher modes can be 
considered as a series of modes associated with vibration of the vertical 
side wall of the diaphragm. The first side wall mode (second tank mode) is 
the dominant mode at low excitation levels. 

A number of conclusions can be drawn from these observations. 

(1) The dynamic behavior of the sloshing tank is highly nonlinear at 
acceleration levels approaching those in flight. 

(2) The slosh modes of the tested tank appear to be highly associated 
with the shape of the diaphragm. The implies that the modes 
should be a strong function of any parameter that reflects a 
change in diaphragm shape. Examples are ullage volume and the 
vertical acceleration field, which effects the hydrostatic fluid 
pressure on the sides of the protruding diaphragm at low ullages. 

(3) Because the primary mode at high acceleration levels is asso- 
ciated with a different, deformed diaphragm shape, one can con- 
clude that the measured higher order modes (which correspond 
to the static undeformed shape) can not occur simultaneously 
with the primary mode. In fact, there probably exists a different 
set of higher order modes associated with the deformed diaphragm 
shape. 

(4) Because the high-level response was significantly different from 
the low-level response, future testing needs to be conducted at 
flight levels to insure that data are applicable to flight analyses. 

(5) During the test program the prototype tank was subjected to a 
number of hours of excitation at the 0. 5 -_g_ level. This level 
slightly exceeds anticipated qualification test levels in the 
critical frequency range below 10 Hz. Because the tank and 
diaphragm exhibited no failures during the test program, an 
important additional conclusion is that the prototype diaphragm 
is structurally quite satisfactory. 
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VI. DYNAMIC PROPERTY DETERMINATION 


A. Parameter Estimation From Nyquist Plots 

Though the quadrature response plots presented in Section V graphically 
display the resonant frequencies, there are other data formats that more 
readily allow the extraction of slosh mass and damping information. The 
Nyquist phase plane plot has been found to provide an ideal format for com- 
plete dynamic property determination and was adopted in this study. Nyquist 
plots of the slosh response at each of the excitation levels are provided in 
Figs. 8 through 11. The ordinate of these plots is the same normalized 
quadrature (Quad) response that was plotted in Fig. 6. The abscissa is the 
normalized coincident (Co) response, in other words, the component of the 
fluid reaction force that is in phase with the driving acceleration. Precisely, 
it is the in-phase component of the fluid reaction force (total in-phase force 
minus rigid fixture and tank accelerating force) divided by the total fluid 
mass times the acceleration level. The Nyquist plots thus can be interpreted 
as phase plane plots of the reaction force amplification factor, with frequency 
as the parameter. 

The Nyquist plot of an ideal linear- single-degr ee-of-fr eedom system 
with viscous damping is a generally circular plot as shown in Fig. 12. The 
plot originates at the point (1,0) at zero frequency and circles to the origin 
as the frequency increases indefinitely. As noted previously, the natural 
frequency corresponds to the frequency near the peak quadrature response. 
The damping is related to the rate of change of frequency along the plot near 
the peak quadrature response, and with the damping known, the system mass 
is determined by the diameter of the circular plot. 

For multiple-degr ee-of-fr eedom systems the Nyquist plot is the sum of 
the circular plots corresponding to each mode, and it generally consists of a 
number of loops as indicated in Fig. 13. With a fluid system each loop cor- 
responds to a slosh mode and the entire plot is shifted to the right by the 
fraction of the total fluid mass (the residual mass) that does not take part in 
the vibration. If the individual modes are reasonably separated so that the 
loop for each mode is well defined, then the dynamic parameters for each 
mode can be estimated from the characteristics of its loop. The parameter 
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values for the i-th mode are defined below in terms of the symbols defined 
in Fig. 13. 


Natural frequency f 


Damping ratio (y 



Mass fraction (p) = 2yD 


( 1 ) 


The point defined by the frequency f is any point in the neighborhood of the 
local peak quadrature response, and the diameter D is for the circle that 
is the best fit in the neighborhood of the local peak quadrature response. To 
avoid the problem of metric mass units or U. S. customary weight units, the 
fluid mass associated with the i-th slosh mode is defined as a fraction of the 
total fluid mass. The mass fraction (p.) for the i-th mode is thus the ratio 
of the i-th slosh mass to the total fluid mass. 

B. Parameter Optimization by Least Squares 

When the individual modes are closely spaced in frequency, highly 
damped, or otherwise coupled, the loops of the Nyquist plot may not be well 
defined. Under these circumstances, parameter estimation based on the 
relationships in Eq. (1) becomes very approximate. This was found to be 
the case with the Nyquist plots obtained in the current slosh test and pre- 
sented in Figs. 8 through 11. 

To alleviate the problem, a nonlinear least squares curve fitting 
approach to parameter determination was developed. In this approach the 
theoretical Co-Quad response of an arbitrary N degr ee-of-fr eedom system 
was derived as a function of the system natural frequencies, damping ratios, 
and mass fractions (see the Appendix). At a specific driving frequency f 
the theoretical response is given by the following expressions: 


Co 




wff 1 - °t) 
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and 


where 
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-2p. £2. Y- 
M 11 


(i . ^ 2 ) + 40 f Yi 2 


F/x = 



n. = 

i 


m. = 
1 

m T = 
f = 

f. = 

l 


driving force/driving acceleration 
mass fraction of i-th mode = m./rrm 
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total fluid mass 
driving frequency 
natural frequency of i-th mode 


U) 


Using the least squares approach, an error function was defined by summing 
the square of the difference between the theoretical and measured response 
over the set of p driving frequencies used in the test. More specifically, 
the error function is a function of the theoretical system state vector (a) and 
is given by: 
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The state vector a, which minimizes Eq. (3), defines the best fit dynamic 
parameter values. 

To allow its application to future programs in addition to the current 
study, the above approach was implemented in the form of a general purpose 
computer program utilizing the variable metric method (Ref. 3) for the func- 
tion minimization. The program NYQFIT was coded with a capacity of 1 to 
10 modes (3 to 30 parameters) and up to 100 frequency response points. In 
addition, either structural or viscous damping was allowed. Computer costs 
were in the one -to-two-dollar range for cases with 5 degrees of freedom and 
around 30 frequency response points. 

Experience with the program indicates that it is a very effective tool 
for Nyquist plot system analysis. One of its primary advantages is in indi- 
cating the response of the best ideal linear system approximation to the mea- 
sured system. The difference between the closest theoretical system 
response and measured response is useful in indicating missing modes, sys- 
tem nonlinearities, and possible data errors. The dynamic parameters 
(natural frequencies, dampings, and mass fractions) associated with the 
closest theoretical system are, in themselves, useful in quantifying the 
measured system’s characteristics. 

C. Parameter Estimates for Tested Tank 

The relationships defined in Eqs. (1) were used to provide initial esti- 
mates of the system parameters from the measured Nyquist plots. The 
program NYQFIT was then used to optimize the parameters to obtain the 
best fit between the theoretical response and the measured response. Nearly 
100 runs were made assuming both viscous and structural damping, and 
using various numbers of modes and various frequency ranges. Most work 
concentrated on the 0. 1-, 0. 3-, and 0. 5-g_ excitation level data. 

Figures 14 and 15 represent initial Nyquist plot fits for the 0. 05 and 
0. 1 -j* data. The parameter values corresponding to the computer generated 
curves are indicated in the figures. The lack of agreement between the best 
fit theoretical response and the measured response indicates considerable 
system nonlinearity in the first two modes. In the measured response the 
first mode resonance appears to partially exclude the second mode resonance. 
Note that the loop associated with the measured first mode response is not 
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pulled as tight by the increasing second mode response as it is in the theo- 
retical response. 

In an attempt to improve the definition of the 0. 1 -g level first mode 
parameters, the first mode only was fitted with various single-degr ee-of- 
freedom models. This work indicated that the first mode natural frequency 
was about 3. 2 Hz, and that the mass fraction was about 15%. When the data 
were incorporated into the overall 4 mode model, the 15% mass fraction was 
found to be inappropriate. As a final iteration the first mode natural fre- 
quency and mass fraction were held fixed while the remaining 10 parameters 
in the 4 mode model were optimized. This was done using various values 
for the first mode mass fraction. Typical results for the constrained 
optimization are presented in Fig. 16 for a first mode mass fraction of 25%. 
Together the data in Figs. 15 and 16 provide a quantitative estimate of the 
system parameter values, and the degree of system nonlinearity at the 0. 1 -j* 
excitation level. 

Attempts to fit the 0. 3 -g excitation level data with a single linear 
dynamic model were much less fruitful. It was quickly established that a 
single multiple-degree-of-freedom model could not be assumed across the 
entire frequency range. As described in Subsection V-B the diaphragm in 
the vibrating tank developed a completely different configuration during the 
primary 4. 8-Hz resonance than it had outside the resonance. Accordingly, 
different models had to be used in the two frequency regions. 

Figures 17 and 18 describe the results of fitting one- and two-degree- 
of-freedom models, respectively, to the primary 4. 8-Hz resonance. As 
seen in Fig. 17 the single-degr ee-of-fr eedom model was unable to account 
for the response pattern in the region around 3 Hz. Based on the 0. 1 -£ level 
response, which indicated a subdued resonance at about 3 Hz, it was con- 
cluded that the 0. 3-_g_ level response probably contained the same resonance 
hidden by the dominant 4. 8-Hz mode. As indicated in Fig. 18 the two- 
degr ee-of-fr eedom model provided an excellent match for the measured 
data. However, the difference between various two-degr ee-of-fr eedom fits 
indicated that the subdued mode is not well defined. This is probably due to 
the fact that the diaphragm is changing shape in this region as the primary 
mode develops. 
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As described in Subsection V-B, the diaphragm jumps back to its 
original dimpled equilibrium shape as the frequency approaches 5 Hz. Fig- 
ure 18 also contains a typical result of fitting 4-degree-of-fr eedom models 
to the response above 5 Hz. The first of the four degrees of freedom was 
used to account for the state of the tank at 5 Hz and was varied parametrically. 
Figure 18 represents the case where the first mode natural frequency and 
mass fraction were constrained to their best estimate values. The program 
was allowed to select the first mode damping as required to reach the start- 
ing point of the 7-Hz mode on the Nyquist plot. 

Considering the 7-Hz mode, it appears to correspond to the 9-Hz mode 
at the 0. 1 - £ level, reduced in frequency because of nonlinearities with exci- 
tation level. The high-order modes are similarly reduced in frequency and 
are generally quite similar to those indicated for the 0. 1-g excitation level 
in Fig. 1 6. 

At the 0. 5 -£ excitation level the primary 4. 8-Hz mode that developed 
at the 0. 3-£ level is essentially unchanged. Figure 19 represents the best 
fit to the 0. 5-g data. Except for the presence of the small resonance around 
9. 6 Hz, which is considered to have been a fixture resonance, the response 
is quite similar to the response at the other levels. The primary difference 
is increased damping and slightly lower slosh frequencies. 

Figures 20, 21, and 22 summarize the approximate variation of the 
modal parameters with increasing excitation level. 

VII. SUMMARY AND RECOMMENDATIONS 

During the test program a 23-cm- (9-in. -) diam mercury propellant 
tank was tested with a 17. 5% ullage volume. Emphasis was placed on estab- 
lishing the general slosh behavior as a function of excitation level, and on 
deriving a computerized data reduction method. Natural frequencies, damp- 
ings, and modal masses were determined for the lateral slosh modes at 
sinusoidal excitation levels between 0. 05 and 0. 5 £ peak. These data were 
reduced with the aid of a Nyquist plot least squares parameter fitting com- 
puter program that was coded as a general purpose tool suitable for use in 
future test programs. 
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The slosh response was found to be very dependent on the stiffness and 
deformed shape of the tank*s positive expulsion diaphragm. As a result the 
response differed significantly from the theoretical response for the same 
fluid volume with a free surface. The response was also found to be very 
nonlinear with excitation level. At excitation levels above 0. 2 _g, the tank 
developed a new primary resonance that corresponded to a new deformed 
diaphragm shape different from the low-level equilibrium shape. Above this 
natural frequency the deformed shape jumped back to the low-level equilibrium 
shape. As a result the higher modes were less affected by the excitation 
level. 

Because the current work concentrated on a single test configuration, 
the understanding of the slosh dynamics is not complete. Further testing 
will be required in the future to identify the significance and understand the 
effects of varying other slosh parameters. At this time the effect of tank 
size, ullage volume, ullage fluid, ullage pressure level, bladder stiffness, 
static acceleration field, and excitation direction is not adequately understood. 

Some insight into the effect of larger ullages and more flexible dia- 
phragms is given by Stofan and Pavli in Ref. 4. Their tests were conducted 
with a similar 24-cm-diam tank with diaphragm thicknesses of 0. 25 and 
0. 5 mm, and ullage volumes of 50% and around 80%. As might be expected 
their data indicate much closer agreement with free surface slosh theory 
than the data presented here for a 1. 6-mm thick diaphragm. However their 
data were obtained at constant excitation levels of 1. 25 mm (0. 05 in. ) and 
2. 5 mm (0. 10 in. ). This corresponds to excitation levels of only about 0. 01 
and 0. 02 g at 4 Hz. As indicated by this study, the results could be signifi- 
cantly different at flight levels that approach 1 j*. 
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APPENDIX 


DERIVATION OF THEORETICAL CO-QUAD RESPONSE 

As indicated in Fig. A-l the tank/fluid system is considered equivalent 
to a set of uncoupled single-degree-of-fr eedom oscillators plus two rigid 
masses. Each oscillator represents a slosh mode and is defined by its stiff- 
ness, mass, and damping. The rigid masses represent the mass of the tank 
and its moving support, and the residual mass of the fluid, which does not 
take part in the sloshing. Though the rotary inertia of the residual mass and 
the equivalent vertical positions of the slosh masses, noted in Fig. A-l, are 
of interest in a complete model of the fluid, only the lateral slosh behavior, 
which is independent of these parameters, is considered. For the total 
response of the model with both pitch and lateral excitation, including the 
height and inertia parameters, see Ref. 5. The following nomenclature is 
adopted: 

m = viscous damping factor for i-th mode 
F(t) = total lateral excitation force 

= scaler amplitude of sinusoidal excitation force 
F = scaler amplitude of sinusoidal fluid reaction force 
g. = structural damping factor for i-th mode 

i = index denoting the i-th slosh mode 

j = J-L 

k. = stiffness of i-th mode 

i 

m. = modal mass of i-th mode 

i 

m^ = total fluid mass 

iyIq = residual fluid mass = m^ - 2 rm 
M = mass of tank and related fixturing 
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n = number of significant slosh modes 
t = time 

x. = lateral displacement of mu relative to the tank 
Xq = absolute displacement of tank relative to ground 
= scaler amplitude of sinusoidal x^ 

X = scaler amplitude of sinusoidal x^ 

X = scaler amplitude of sinusoidal tank acceleration x^ 

•y. = damping ratio for i-th mode = (c/c ) 

1 i 

u. = mass fraction for i-th mode = m./rrm 
i i T 

co = frequency of sinusoidal excitation 

ecu = natural frequency of i-th mode 

O,, = frequency ratio of i-th mode = oo/uu 

X = summation over i for i = l,n 
i 


Writing Newton’s second law for the lateral motion of the total tank/ 
fluid system gives 


F(t) = (M 


+ m o ,x o + Xl m i (i o +i i : 


(A- 1 ) 


and for the i-th slosh mass gives 


-k.x. - c.x. = m.{x„ + x. ) 

li li l 0 l 


or 


m.x. + c.x. + k.x. = -m.x„ 
li li ii lO 


(A-2 ) 
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For steady state sinusoidal tank motion: 


F(t) = F t exp ( j cot ) 

Xq = X exp (jwt) 

x . = X. exp (jwt) 

x . = jwX i exp (jwt) 

2 

x . = -w X. exp (jwt) = X. exp (jwt) 


Substituting into Eq. (A-2) gives 


X. 

l 


2 Y 

co m.X 

l 


-w m. + iwc. -f k. 

l J l l 


Substituting Eqs. (A-3) and (A-4) into Eq. (A-l) and factoring 
acceleration gives: 


F rp exp (jwt) = 


M 4- m.Q + 


E m i + XI 


2 2 
m. w 

i 


-w m. + iwc. + k. 
11 1 y 


X 


or 


T 

— - M = m + 
X 1 
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( A-3) 


(A-4) 


out the tank 


exp (jwt) (A- 5) 


(A-6) 


Now define 


F = F t - MX 


z = jh 

i m. 

1 


jc_\ _ c i 

' i l c ) 2m. to. 

V c /. 11 


Q. 

1 


CO 


CO. 

1 


h = 


m. 

l 


i m r 


Substituting the definitions of Eqs. (A-7) into Eq. (A-6) gives 


F 


X m r 


E M-.a* 

LLJ 

i - + ziv.a 


2jy. 

i J "i i 


Separating Eq. (A-8) into its real (Co) and imaginary (Quad) 
nents gives: 


Co 


. -rp . r — > p.a (i - a \ 

AM = ! + V — EA — iL. 

\xm T ) i (l - E) + 4 vM 


Quad 


mu = y — 

' Xm T / i (l - ^ + 4v^£2? 


(A-7) 


(A-8) 


compo- 


(A-9) 


(A- 1 0 ) 
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which is the desired solution for viscous damping. For structural damping 
the result is identical to the above except that the viscous damping factor y* 
is replaced by 



(A- 1 1 ) 


where g^ is the structural damping factor for the i-th slosh mode. 
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Fig. 1. SEPST Program propellant tank design (propellant-loaded position): 
(a) line drawing; (b) disassembled mercury propellant tank 
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Fig. 2. Test tank mounted in sloshing apparatus 
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Fig. 4. Test tank with water (17. 5% ullage) 
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Fig. 3. Shaker control and data processing console 
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Fig 


. 5. Static diaphragm shapes 
(a) 90% ullage; (b) 80% ullage; 


with mercury for various ullages 
(c) 70% ullage; (d) 60% ullage 
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Fig. 5 (contd) 
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Fig. 5 (contd) 
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Fig, 5 (contd) 
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QUADRATURE RESPONSE, Quad 


iN-PHASE RESPONSE, Co/-!—' 
\ m T> 


- 0.4 - 0.2 0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



Fig. 8. Nyquist plot for 0.0 5 g_ excitation level 



Fig. 9. Nyquist plot for 0. 1 £ excitation level 
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(N-PHASE RESPONSE, Co 


CO 

00 

vO 


CO 

CO 



CO 

CO 


Fig. 13. Typical Nyquist plot for multiple- degree-of-f reedom system 






QUADRATURE RESPONSE, Quad (77— 


IN-PHASE RESPONSE, Co/-^- N 
v m T/ 



Fig. 14. Initial linear system fit for 0. 05 £ data 


IN-PHASE RESPONSE, Co 

\ m T / 



Fig. 15. Initial linear system fit for 0. 1 g_ data 
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QUADRATURE RESPONSE, Quad (77— 


IN-PHASE RESPONSE, Co/77^- 
V* m T, 


0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



Fig. 16. 


Optimized linear system fit for 0. 1 g data 


IN-PHASE RESPONSE, Co/-!— 
\ m T/ 


0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



Fig. 17. Initial single-degree-of-freedom fit for 0. 3 £ primary resonance 
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QUADRATURE RESPONSE, Quad ( T f— \ QUADRATURE RESPONSE, Quad 



Fig. 18. Optimized linear system fit for 0. 3 data 


IN-PHASE RESPONSE, Co/- 1 --] 

\ Xm T/ 

0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



Fig. 19. 


Optimized linear system fit for 0. 5 £ data 
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Fig. 22. Mass fraction versus excitation level 
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Fig. A-l. Dynamic model of tank/fluid system 
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